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ABSTRACT 

We analyse whether a stellar atmosphere model computed with the code CMFGEN 
provides an optimal description of the stellar observations of WR 136 and simultane¬ 
ously reproduces the nebular observations of NGC 6888, such as the ionization degree, 
which is modelled with the pyCloudy code. All the observational material available (far 
and near UV and optical spectra) were used to constrain such models. We found that 
even when the stellar luminosity and the mass-loss rate were well constrained, the 
stellar temperature T* at r = 20, can be in a range between 70 000 and 110 000 K. 
When using the nebula as an additional restriction we found that the stellar models 
with r* ~ 70 000 K represent the best solution for both, the star and the nebula. 
Results from the photoionization model show that if we consider a chemically homo¬ 
geneous nebula, the observed N“'"/0^ ratios found in different nebular zones can be 
reproduced, therefore it is not necessary to assume a chemical inhomogeneous nebula. 
Our work shows the importance of calculating coherent models including stellar and 
nebular constraints. This allowed us to determine, in a consistent way, all the physi¬ 
cal parameters of both the star and its associated nebula. The chemical abundances 
derived are 12 -|- log(N/H) = 9.95, 12 -|- log(C/H) = 7.84 and 12 -|- log(0/H) = 8.76 
for the star and 12 -|- log(N/H) = 8.40, 12 -|- log(C/H) = 8.86 and 12 -|- log(0/H) = 
8.20. Thus the star and the nebula are largely N- and C- enriched and 0-depleted. 
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1 INTRODUCTION 

NGC 6888 is an emission nebula associated with the Wolf- 
Rayet star WR 136 with spectral type WN6(h), and it is a 
typical example of a ring nebula resulting from the continu¬ 
ous interaction between a strong stellar wind and the inter¬ 
stellar medium. Projected on the sky, this nebula seems to 
be a filamentary and clumpy gaseous ellipse emitting mainly 
in Ha and optical [N II] lines, with an axis-size of 12 x 18 ar- 
cmin ( Chu, Treffers, fc Kwitte'r||1983 1, and it is surrounded 
by a homogeneous sphere which emits mainly in [O III] 5007 


A ( Gruendl et al.|2000 Moore, Hester fc Scowen|2000 1. The 


combination of its filamentary morphology and its nitro¬ 
gen enhanced abundance suggests that the ejections, coming 
from a previous evolutionary stage of the star WR 136, have 
contributed in an important way to the formation of the ring 


nebula that we see today ( 

Esteban & Vflchez|1992, hereafter 

EV92||Fernandez-Martm et al.|2012, hereafter FM12||Mesa- 

Delgado et al.|2014, hereafter M-D14 

1. 


Both the nebula NGC 6888 and its central star WR 136 
have been already studied in several wavelength ranges from 
X-rays to radio emission. For example, the stellar parame¬ 
ters, such as the stellar temperature, mass-loss rate and ter¬ 
minal velocity, have been calculated by|Crowther fc Smith 


(19961 and Hamann, Grafener & Liermann (2006, hereafter 


HGL061 by using stellar atmosphere codes and comparing 


the results with observed ultraviolet, optical and infrared 
stellar spectra; while the wind velocity law and dumpi¬ 


ness were studied by Ignace, Quigley & Cassinelli (20031. 
[Skinner et al.| ( |2010| ) showed that WR 136 emits in X-ray 
energies. In particular, these authors found soft emission 
(kT ~ 0.56 keV) coming from outside the wind and hard 
emission {kT ~ 2.64 keV) coming from the gas inside the 
wind. 

The complete hot bubble in the interior of NGC 6888 
was observed with the Space Telescope Suzaku (jZhekov'^ 


Park 20111, confirming the existence of nitrogen enriched 


E-mail: jperez@astro.unam.mx (JRP) 


X-ray emitting plasma which consists of two thermal com¬ 
ponents, one emiting soft X-rays and the other one, hard 
X-rays. Based on this, the authors proposed a new X-ray 
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emission mechanism based on the thermal interaction be¬ 
tween heavy particles, given that the known mechanisms 
are not capable to explain simnltaneonsly these two X-ray 
emitting components. However, as these authors explain, it 
is necessary to have better observational restrictions to con¬ 
firm this idea. 

On the other hand, complete hydrodynamical simnla- 
tions performed by Toala & Arthur (20111, attempting to 
reproduce the X-ray observations, indicate that the nebular 
morphology plays a key role in the presence of X-rays in ring 
nebulae. These authors conclude that the X-ray emission is 
higher in bubbles with shell-like morphology than the one 
present in clumpy filamentary bubbles. However, as the au¬ 
thors indicate, the stellar evolutionary models are not able 
to reproduce simultaneously the observed stellar parameters 
and the spectral signatures of the observed X-ray luminosity. 

The physical properties in different zones of NGC 6888, 
such as electronic temperature and density, expansion ve¬ 
locity and chemical composition, have been investigated by 
several authors since [Kwitter] ( |1981[ ) to recently [Stock k\ 
Barlow| 120141. EV92 made a chemodynamical study and 


they found that NGC 6888 has a chemical composition re¬ 
vealing the material processed in the interior of a star with 
an initial mass between 25 M© and 40 M©. These results 
were confirmed by M-D14 who detected for the first time 
the C II 4267 A recombination line, allowing to estimate the 
C abundance. In this way they presented a complete study 
of the trace of CNO cycle in the nebula. 

FM12 made a very complete study of NGC 6888 
through integral field spectroscopy observations performing 
ID and 2D analyses. Their results show indeed that NGC 
6888 has a complex structure with shells and filaments mov¬ 
ing with different velocities. In particular they proposed an 
“onion-like” morphology composed mainly by three shells. 
Concerning the plasma properties, they found that the elec¬ 
tron density ranges from less than 100 to 360 cm“® and 
the electron temperature presents variations from ~7 700 to 
~10 200 K. In addition, they reported important variations 
in the nitrogen abundance across the nebula, the N/O ra¬ 
tio being much higher in the inner region. Recently jStockj 
jfc Barlow ( 2014| ) derived N/O values from the infrared [N 
III] 57-/rm and [O IIlj-88 /rm lines and their results would 
indicate that the edge region of NGC 6888, over the major 
axis, has lower N abundance with respect to the inner region 
near the star. 


All the previous works focus only on the modelling of 
the star or the nebula but not both simultaneously. Fur¬ 
ther, when a single study of these systems is made (star or 
nebula) it is difficult to assign a value to the physical param¬ 
eters only based on generic data. This is why we focus in a 
semiempirical method based in the self-consistent modelling 
of the star and its associated nebula. 


It is logical to proceed in this way because the nebula 
associated to the star is naturally an additional restriction 
to the stellar atmosphere model, because the nebula can 
tell us about the ionizing ultraviolet photons coming from 
the star and then reprocessed and reemitted with low energy 
through recombination and collisional processes. In this way 
we can realize how well a stellar model can produce the ion¬ 
ization stage of the nebula and for this purpose we build a 
photoionization model of NGC 6888 based on spectroscopic 



Figure 1. Slits positions where the observations were performed 
on NGC 6888 including the positions of M-D14 labeled as M-D, 
B2 by FM12 and the position of EV92. Red filter image from the 
STScI-DSS I/II server. 


observations where the modelled stellar spectral energy dis¬ 
tribution (SED) is used as the ionization source. 

This methodology will give us the advantage to obtain 
more constrained and coherent physical values for all the 
components of the entire system than those resulting from 
a single study of each component as it was made until now. 
The results from this work will help to better constraint 
other kind of studies as models of stellar evolution, astro- 
physical bubble formation and, in particular, the presence 
of X-rays in the bubbles since this hot gaseous component 
imposes an additional and important restriction to under¬ 
stand the whole system, given that is there where the stellar 
wind energy is stored. 

This paper is structured as follows: In Section 2, ob¬ 
servations and data reduction process are described. Sect. 3 
contains the information of the line flux measurements and 
the error estimates. In Sect. 4 we describe the stellar and 
nebular models as well as our methodology to constraint the 
free parameters. A discussion about our best self-consistent 
model is presented in Sect. 5. Finally, in Sect. 6, we present 
our conclusions. 


2 OBSERVATIONS AND DATA REDUCTION 


The main analysis in this work is based on the optical ob¬ 
servations carried out the night 1997 August 28, for the 
star WR 136; and the 2011 August 31 and September 1, 
for the nebula and the standard star HR 7596. We used 


the REOSG echelle spectrograph (Levine & Chakrabarty 
|1994p attached to the 2.1-m telescope at the Observatorio 
Astronomico Nacional (OAN) in San Pedro Martir, B. G. 
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Table 1. Log of observations 


Observation ID 

A. R 
J2000 

Dec. 

J2000 

Date 

Exp. Time 

(s) 

Slit dimensions 
(arcsec) 

NGC 6888A 

20:12:36.5 

38:26:56.5 

1st sept. 2011 

4x1800 

4x13.3 

NGC 6888B 

20:12:03.7 

38:25:07.8 

30 Aug. 2011 

4x1800 

4x13.3 

NGC 6888C 

20:11:51.2 

38:21:45.8 

30 Aug. 2011 

4x1800 

4x13.3 

WR 136 

20:12:06.5 

38:21:17.7 

28 Aug. 1997 

3x300 

9.3x13.3 

HR 7596 

19:54:44.7 

00:06:25.1 

1st sept. 2011 

3x420 

9.3x13.3 


Mexico. With this spectrograph we obtained high-resolution 
spectra [R ~ 20000) with a wavelength coverage from 3700 
A to 7300 A. 

Three different positions of the nebula were observed 
using a slit with size 4 x 13.3 arcsec oriented east-west. For 
each position we obtained four spectra with an exposure 
time of 1800 s each. The slit positions are shown in Figurej^ 
Together with slit positions of other authors that will be 
described in the following sections. 

To perform the spectroscopy of WR 136 and the stan¬ 
dard star the slit size was 9.3 x 13.3 arcsec. For each star 
three spectra were obtained, with a exposure time of 420 s 
in the case of the WR star and 300 s for the standard star. 
The log of observations is shown in Table 

After the bias correction was made, the corresponding 
different exposures were combined and then background and 
cosmic rays were subtracted using the filter/cosmic task of 
the astronomical software MIDAlQ Then, once the wave¬ 
length calibration was performed (using a ThAr lamp), the 
orders of the echelle images were extracted. 

To correct the extracted orders by the blaze function 
and to make the flux calibration in the case of the nebular 
spectra, we modelled the star HR 7596. First the individual 
orders were normalized to a curve which is theoretically the 
convolution between the blaze function and the sensitivity 
function. In this way, the normalized orders were merged to 
obtain a complete normalized stellar spectrum to be mod¬ 
elled with the code SPECTRUM ( [Gray] 199^ . This code as¬ 
sumes a LTE atmosphere with a plane-parallel geometry and 
is adapted for modelling standard stars with a spectral type 
between B and M. The spectrum can be calculated in a nor¬ 
malized way or the model can give the emergent flux per 
wavelength unit. First we obtained the normalized modelled 
spectrum to directly compare it with the observed spectrum 
and once this one was reproduced, the not normalized spec¬ 
trum was calculated. 

The model spectrum of the standard star was then red¬ 
dened considering the interstellar and atmospheric extinc¬ 
tion as well as the radiation dilution. In this way, we con¬ 
structed the synthetic standard star spectrum. Then, the 
extracted orders of the standard star observation can be 
divided by this spectrum to obtain the single order func¬ 
tion where the information about the instrumental sensi¬ 
tivity and flux calibration is contained. To correct for the 
CCD sensitivity and to flux-calibrate the nebular spectra, 
these orders were divided by their corresponding single or¬ 
der function. Finally the corrected extracted orders were 


^ http://www.eso.org/sci/data-processing/software/esomidas 


analysed as individual spectra. These resulted spectra have 
a total exposure time of two hours. 

In the case of the orders from the WR 136 frames, we 
proceeded with the rest of the reduction similarly as we 
made with the standard star and given that observations of 
the standard star and WR 136 were not made at the same 
date, the resulting optical spectrum was just normalized to 
the continuum. 

In addition to the optical spectra of WR 136, we an¬ 
alyzed a spectrum in the UV region between 900 A and 
1200 A obtained from the FUSE satellite observational pro¬ 
gram C097 where some interesting and useful lines, such 
as S V A1063 and P V All 18 lie. To extend the wave¬ 
length range to 3300 A we retrieved several reduced lUE 
spectr^ To improve the signal to noise ratio we com¬ 
bined the individual lUE short-wavelength high-dispersion 
spectra SWP57618, SWP57622, SWP57625, SWP57626, 
SWP57630, SWP57635, SWP57638, SWP57645, SWP57647 
and SWP57663; and the individual lUE long wavelength 
high dispersion spectra LWP32398, LWP32407, LWP32409, 
LWP32413, LWP32416, LWP32418, LWP32423, LWP32440, 
LWP32441 and LWP32446. 


3 FLUX MEASUREMENT AND ERROR 

ESTIMATES 

In Figure we show the resulting Hq and [N II] 6548, 6584 
A lines present in each spectrum of NGC 6888. In this figure 
it is apparent that the spectra show line emissions coming 
from three components at different velocities. This is more 
obvious for the zone B in whose spectrum we can clearly dis¬ 
tinguish these components located at around —45, 15 and 65 
km s“^. According to the wavelength calibration performed, 
the radial velocities have an accuracy of ±12 km s“^. 

In the spectrum of position C only two components are 
present with radial velocity of —57 and 12 km s“^. Probably 
the redshifted one is also present but very weak. In the case 
of position A, corresponding to the nebular edge, it is only 
evident the more intense blueshifted component at around 
—36 km s“^. We can see that the line profiles are not gaus- 
sian at all, they exhibit a wing on the red side where surely 
the others components, not resolved, are contained. 

The radial velocities of the blueshifted and redshifted 
components agree with the values in the literature. For ex¬ 
ample, in a position near the star, EV92 detected two com¬ 
ponents at —64 and 78 km s”'^ and M-D14 reported two 

^ see https://archive.stsci.edu/iue/file_formats.html for more in¬ 
formation 









4 J. Reyes-Perez, C. Morisset, M. Pena and A. Mesa-Delgado 


[N II] A6548 Ha [N II] A6584 









1.0 
0.8 
0.6 
0.4 

0.2 
0.0 

-100-50 0 50 100 

V km s"' 

LSR 




Figure 2. The Ho and [N ii] 6548, 6583 lines are shown in radial velocity space for positions A (top), B (middle) and C (bottom). The 
velocities are measured respect to the Local Standard of Rest. Spectra at position A only present one component at -36 km s“^, while 
the spectra at position B show two components corresponding to the expanding shell (bluest and reddest components) at -45 and 65 km 
s“^ and a central component corresponding to the H ii region around NGC 6888 at 15 km s“^. Finally in the spectra at position C two 
components are detected at -57 and 12 km s”'^. 


blueshifted components at —60 and —25 km s“^, all of them 
associated with the emission of NGC 6888. Recently, per¬ 
forming a 2D study in a region close to our position A, FM12 
found that the emission of the [N II] 6584 A line comes from 
plasma with radial velocities between —40 to 30 km s“^, in 
agreement with our value at position A. 

The kinematical component at ~15 km s“^ is surely 
associated with the extended H II region around NGC 6888 
and it has been also identihed by EV92 and M-D14 with a 
velocity of 18 km s“^ and 12 km s“^, respectively. 

In general, we found radial velocities which follow the 
trend previously reported, i.e, the highest values appear near 
the star and they decrease as moving away from the centre 
of the nebula. 

For the spectra at positions B and C we measured the 
line fluxes by htting a gaussian prohle to each kinemati¬ 
cal component, but for the spectrum at position A, all the 
emission was integrated for each identified line. The errors 
assigned to our measurements correspond to a two sigma 


value over the mean noise measured at each side of the lines 
detected. 

The fluxes were then dereddened with a logarithmic ex¬ 
tinction coefficient cnp determined from the Ha/ H/3 ratio 
using [Cardelli, Clayton fc Mathis| ( |1989[ ) extinction law. The 
uncertainties in cnp are the result of the flux-error propa¬ 
gation. 

The crr /3 value at position A (0.70) agrees with the mean 
value of 0.76 obtained by FM12 in a region close to this 
observation. In the case of position B the extinction values 
are 0.36, 0.38 and 0.77 for the blueshifted, foreground and 
redshifted components, respectively. The first and second 
values are similar to the value 0.34 obtained by M-D14 for 
a region near our position B. 

The high ch^ value of the redshifted component seems 
to indicate that the emission coming from the back shell 
of NGC 6888 is more reddened surely by the presence of 
inner extinction. This is supported by the extinction value 
for WR 136, Av = 1.73 mag, reported by |van der Huch^ 
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Table 2. Fluxes (normalized to = 100) of the identified lines in the observed spectra, for positions A, B (blue, central and red) 

and C(blue and central). The flux Yup is in units of 10“^^erg cm“^ s"'^. 


A (A) 

Ion 

A 


Be 

Br 

Cb 

Ce 

3726.03 

[O II] 

52.6±14.3 






3728.81 

[OII] 

62.6±9.5 






3970.07 

H, 

16.8±6.1 






4101.73 

As 

28.1±13.0 

19.1ill.O 



24.2i9.4 


4340.46 


47.1±5.7 

48.8il3.4 

26i20.5 

36.9i21.4 

47.7i 3.5 

43i24.0 

4471.47 

He I 

6.Oil.2 

4.2il.7 

3.5i2.4 

2.9i3.3 

7.5i2.1 


4958.91 

[O III] 

23.8i2.4 

97.0i6.2 

72.9i24.1 

101.2i22.8 

52.1i4.5 

77.8i47.5 

5006.84 

[O III] 

67.9i5.6 

289.2il6.3 

210.9i66.4 

250.2i57.2 

154.2i9.2 

230.2il01.5 

5015.67 

He I 

4.7il.8 

e.Oil.S 


44.7i27.2 

4.7i3.1 

6.8i5.3 

5754.64 

[N H] 

3.8il.2 

2.3il.3 



2.3i0.9 


5875.61 

He I 

19.9i2.3 

21.3i2.5 

14.5i6.7 


23.5i2.1 

18ill.9 

6312.10 

[S III] 


2.2il.3 



1.4i0.6 


6548.04 

[N H] 

168.4ill.5 

37.1i2.8 

28.1il0.4 

81.4il8.8 

70.9i3.5 

40.7il7.5 

6562.80 

He 

285.9il8.3 

286.4il4.1 

285.9i84.4 

286.2i58.8 

286il0.8 

286i92.14 

6583.46 

[N H] 

502.1i31.9 

106.8i6.1 

75.5i24.4 

236.6i49.7 

227.0i9.7 

123.6i46.9 

6678.15 

He I 

6.4il.2 

6.4i2.8 

3.2i2.6 

8i5.7 

5.2i0.7 


6716.44 

[S II] 

22.9i2.5 

7.0i2.1 

24.9il0.0 

8.4i5.5 

7.4i0.8 

30.1il3.1 

6730.81 

[S II] 

19.9il.9 

5.7i2.4 

17.8i7.1 

7.1i4.5 

6.3il.2 

19.9i9.6 

7065.17 

He I 

3.2il.O 

6.0i2.6 



6.4i0.8 


7135.80 

(Ar III] 

9.2il.l 

15.6i2.7 

9.0i5.8 

12.4i6.7 

10.Oil.1 



FjT/3 

4.70i0.05 

0.110i0.004 

0.077i0.004 

0.037i0.005 

0.565i0.023 

0.076i0.022 


chp 

0.70i0.02 

0.38i0.07 

0.36i0.11 

0.77i0.17 

0.0 

0.03i0.45 


(20011. Using Rv = 3.1 and Eb-v = 0.692 x cup, as derived 
from the extinction law mentioned above, we find ch /3 = 
0.80 for the star. Given that WR 136 is embedded in NGC 
6888, inner extinction could be contributing to redden its 
radiation. 

Finally, observations of position C do not seem to be 
reddened. The only observations to compare with a neigh¬ 
bour zone are those performed by EV92. They found cnp = 
0.28, which is the lowest value reported in the works cited 
above, indicating this region of NGC 6888 does not present 
high extinction. 

In Table[^the whole set of identified lines with their re¬ 
spective dereddened flux-integrated values are shown. In this 
table we present one set of values for position A, three (B;,, 
Be and Br) for position B, and two (Cf, and Cc) for position 
C, where the subscripts “b”, “c” and “r” correspond to the 
blueshifted, central and redshifted components, respectively. 


4 THE MODELS 
4.1 The stellar model 


To model the observed stellar spectra of WR 136, in the 


UV and optical ranges, the CMFGEN code | Hillier & Miller 


1998) was used. CMFGEN (CoMoving Frame GENeral) is a 


non-LTE line-blanketed atmospheric code designed for spec¬ 
tral analysis of stars with stellar winds. It solves the radia¬ 
tive transfer and statistical equilibrium equations consider¬ 
ing spherical symmetry. This code is adequate to model the 
atmosphere of massive stars with an extended wind that is, 
in general, larger than the hydrostatic radius of the star. 

To start-up CMFGEN code it is required to set the 
stellar hydrostatic structure through input parameters 
such as the gravity and stellar radius or, in an equivalent 


way, to set a velocity law that approaches in depth to 
the hydrostatic structure. In this study we chose the later 
option and used the velocity law given by: 


U(r) = 


Vo + {Voc-Vo){l-R*/rf 

1 -I- {VoIVcore - l)eXp{{Rt, - r]//le//)’ 


( 1 ) 


as proposed by Hillier et al. (2003). In expression (1) Vc 


is the velocity at the hydrostatic stellar radius R *, Vo is the 
photospheric velocity which controls the match between the 
hydrostatic atmosphere and the wind, Uoo is the terminal 
velocity, he// is the scale height of photosphere, in units 
of Rf , which specifies the density structure at low velocities 
and /? is the so-called velocity law exponent. The values used 
for all our models for Wore, W and he// were 1 km s“^, 100 
km s“^ and 0.002, respectively, as it is suggested for WR 
star^ The /3 values used for the models were 1 or 2, which 
are usual values for this kind of stars (Ignace, Quigley & 
|CassinellH2003p . 

Then, once the code calculates the hydrostatic values, 
it solves the radiative transfer equations in the wind zone 
and the result is the emergent flux (as emitted at 1 kpc 
from the observer) as a function of the wavelength. To fit 
the UV spectra, we used for the modelled spectrum a dis¬ 
tance of 1.45 kpc, based on the HIPPARCOS parallax ( |van| 
Leeuwen|2007 1, and a final extinction value Eb-v of 0.52 (or 


chp = 0.75) which agrees with the value by van der Hucht 
| |2001[ | shown in §3. In this way we fit the stellar continuum 
and the absorption feature around 2200 A. For the optical 
spectrum, the model was analysed using the continuum nor¬ 
malized spectrum. 


® See the GMFGEN Documentation and WR models at kook¬ 
aburra.phyast.pitt.edu/hillier/web/GMFGEN.htm 
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Table 3. Stellar parameters 


Model 

rM 

log L 

[Le] 

log M 
[M0yr-l] 

Voo 

[km s“^] 

T. 

[kK] 

D/3 

[kK] 

d 

fv 

X 

Y 

This work starl 

1.50 

5.45 

-4.63 

1550 

110 

55 

2 

0.1 

0.016 

0.95 

This work star2 

2.10 

5.40 

-4.95 

1550 

90 

48.30 

1 

0.175 

0.068 

0.91 

This work starS 

3.38 

5.40 

-4.95 

1550 

70 

45.91 

1 

0.25 

0.12 

0.86 

This work star4 

3.38 

5.40 

-4.95 

1550 

70 

42.26 

2 

0.25 

0.12 

0.86 

Hamann et aL(2006) 

3.34 

5.40 

-4.95 

1600 

70.8 


1 

0.25 

0.12 

0.86 

Crowther & Smith (1996) 

5.00 

5.30 

-3.91 

1750 

55.7 

28 






^ and T* are defined at r = 20 
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Figure 3. Selected lines from the observed spectra to constrain the models. In each plot the dashed black line is the observed spectrum; 
red, green, yellow and blue lines are respectively the starl, star2, starS and star4 models obtained through the CMFGEN code but 
reddened. The optical spectrum is shown as flux-normalized. 


For the stellar models computed in this study, the in¬ 
put values used to constrain the model atmosphere for WR 
136 were the mass-loss rate M, the terminal wind velocity 
Kxd, the stellar luminosity L, the hydrostatic radius of the 
star Rt (at optical depth r equal to 20) and the chemical 
composition. 


According to the stellar wind theory, it is expected that 
the recombination line luminosity L\ is proportional to the 
mass-loss rate through the relation M ~ Vcoy/Tp, therefore 
the Hq recombination line wa s used to adjust the value of 
M ( Lamers fc Cassinelli[|l999 1. 

On the other hand, when we consider an extended stel- 
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Figure 4. General view of the starl, star2, starS and star4 models compared with the observed spectra (dashed black line). From up 
to down: the FUSE FUV spectrum, the lUE UV spectrum and the continuum-normalized optical spectrum. The colors are as in Fig.|^ 


lar atmosphere, it is not possible to assign a surface temper- 
atnre to the star as the temperature varies with the radius, 
therefore the temperature cannot be a global input value 
in the model. However, it is possible to define a tempera¬ 
ture value T* for r = 20 and such a parameter is set in the 
model through the luminosity and the stellar radius with the 
common expression L = AnaRlTf, where a is the Stefan- 
Boltzmann constant. Therefore, for constructing the model, 
the L and R, values were changed to fit the line intensities 
corresponding to consecutive ionization degrees of the same 
element, such as the He I and He II lines. 


While constructing a model and when the He lines tem¬ 
perature indicators (He II A4686 and He I A5876) pointed 
that we were very close to the best fit, we found that the 
predicted N III A4640 line was weaker than the observed 
line, while the predicted N IV A4058 was more intense than 
the observation, implying the temperature should be smaller 
(see Fig.[^. It has been stated in the literature that the ion N 


III is very dependent on the detail of its atom model | Hillier 
et al. ][2003| . We tried to correct this issue by using a more 


complex N III ion model, i.e., more levels were considered 
when decollapsing a superlevel. However, the problem was 


not solved although the fit of these nitrogen lines improved 
slightly. 

The value of the terminal velocity was constrained 
by fitting the blue edge of the saturated C IV A1550 P- 
Cygni profile. The obtained value (1550 km s“^) fits well 
other P-Cygni profiles such as S IV AA1063,1073 and P V 
AA1118,1128 in the far UV spectrum of the star (see Fig.j^ 
and it is in excellent agreement with values reported in pre¬ 
vious works. [Prinja, Barlow Howarth| ( |1990[ ) used the com¬ 
mon diagnostic line, the saturated UV C IV resonance dou¬ 
blet at 1548, 1550 A and obtained Voo = 1605 km s“^; 


Ee- 


nens & Williams 

(1994 

1 derived Voo = 1660 km s ^ by using 

the He I 10830 .A 

line width; HGL06 measured Voo = 1600 


km s ^ by fitting their model to their observed spectrum; 


Crowther & Smith 


(19961 determined a Voo = 1750 km s 


finally, through an IR study, Ignace et al. (20011 used the 
Ca IV 3.207 fim forbidden line width, reporting Voo = 1490 
km s“^. As our fit was done measuring the velocity at the 
edge of the satured component of the blue wing, instead of 
measuring it where the wing reaches the continuum, this 
renders a terminal velocity ~ 200 km s“'^ lower than when 
measured at the continuum. This excess in velocity can be 
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explained invoking a turbulent velocity or by the clumping 
factor fv of the wind. 

The chemical composition in the models was initially 


set at typical WN composition (van der Hucht, Cassinelli 
fc Williams||1986| and it was modified to generate some of 
the models shown here. The final H and He abundances (X, 
Y by mass-fraction) employed for all the stellar models are 
shown in Table while for the other elements such as C, 
N, O, S and Fe, the abundances derived from CMFGEN are 
shown in Tablej^to compare then with the resulting nebular 
abundances from the modelled nebula. 

As the modelled spectral features do not respond lin¬ 
early to the input parameters, once a value is found it is 
not necessarly fixed but only constrained while changing the 
other input values. Therefore we calculated a grid of models 
in order to find the more appropriated models which repro¬ 
duce most of the observed spectral signatures. 

To properly fit the luminosity through the stellar con¬ 
tinuum without changing the stellar temperature, it was nec¬ 
essary to use the transformed radius 


Rt — R> 


f 14o/2500 
VM/D/10-4 


2/3 


( 2 ) 


(where D = fy^) because it has been shown that two mod¬ 
els with the same Rt and T* values have the same spectral 
behaviour ( [Schmutz, Hamann fc Wessolowski|1989| ). 

From the grid of models calculated by varying the input 
parameters we obtained starl which fits very well most of 
the stellar spectral characteristics. This model is presented 
in Table|3l In addition we also show in this table models from 


the literature, as HGL06 and Crowther &i Smith (19961 ones. 


The mass-loss rate value of starl is slightly higher but of the 
same order as HGL06 value and it is lower than the value 
published by [Crowther fc Smith (19961. We are conhdent in 
the mass-loss rate value of our models because as it can be 
seen in Figure the recombination lines, like the Balmer’s 
series in the optical spectrum, are very well fitted. 

In Figure [^we can see that the He II A1640 and A5411 
lines are overestimated in starl only by a factor ~ 1.3. It 
is important to say that it has been stated in the literature 
that the He lines in an observed spectrum cannot be fitted 
simultaneously ( Hillier et al.|[2003 l. From Figure]^ we can 
also see that starl reproduces in general very well the ob¬ 
served spectra from the UV to the optical wavelength range. 

The main difference between starl and the model re¬ 
ported by HGL06 is the stellar radius i?* which in our case 
has a low value of 1.5 Rq (compared to the value of 3.34 Rq 
by HGL06) leading to a higher stellar temperature T* of 110 
kK. Additionally the chemical composition, represented by 
the values of X and Y, is different. These differences have 
important implications when modelling the nebula around 
the star (see §4.2.2) and are indicating that it could exist 
more than one solution for a stellar model of a WR star 
as both, starl and HGL06 model, reproduce similarly well 
the stellar spectrum (see Fig. 68 by HGL06). This apparent 
degeneracy has already been remarked in a study of WR 


stars modelling by Hillier (19911. In this work the origin of 


this degeneracy is attributed to our misunderstanding of the 
dynamics in the inner wind regions. 

To investigate this degeneracy we searched for other 
CMFGEN models that could fit the stellar observations. The 


results are plotted in Figure where the general behaviour 
of the models is compared with the observed spectra. 

In a simple manner we used interpolated values between 
the model starl and the HGL06 model for the radius, the 
clumping factor and the hydrogen and helium fractions by 
mass, X and Y. The result is model star2 which has an 
effective temperature of 90 000 K (its characteristics are 
presented too in Table|^. Model star2 reproduces the stellar 
spectrum similarly well as starl. It is presented in Fig.|^ in 
green. 

We computed two additional models. The first is star3 
which basically uses HGL06 model parameters but it was 
calculated with the CMFGEN code. Model star3 reproduces 
very well the results of HGL06 model, which shows that 
CMFGEN and the Potsdam WR model atmosphere code are 
equivalent. The last one is model star4 which is almost the 
same as star3 but calculated with the /3 value equal to 2, 
leading to a lower T 2/3 value because the wind region be¬ 
comes more extended. Models star3 and star4 reproduce 
fairly well the observed stellar continuum and many spec¬ 
tral features (lines yellow and blue in Fig. although the 
emission lines indicative of the stellar temperature are worse 
fitted than with starl (see Fig.[^. All the stellar model char¬ 
acteristics are presented in Table 

From this, it is then clear that there is not a unique stel¬ 
lar model reproducing the stellar spectrum. As said above, 
we are in presence of a degeneracy occurring when modelling 
WR stellar atmospheres and evidently more constraints are 
required to decide for the best model. Additional constraints 
can be found by analyzing the ionized nebula. 


4.2 The nebular model 

4-2.1 Abundances and plasma analysis 

The starting point to make a photoionization model for NGG 
6888 is through the spectral analysis of the observations 
taken for different zones on this nebula. 

As we said before, the fluxes obtained from the spec¬ 
trum at position A contains not only the nebular emission 
but also the foreground emission. This foreground contribu¬ 
tion is known from the spectra at positions B and G, where 
the flux values corresponding to the 15 km s“'^ component 
are similar. This enforces the assumption that these compo¬ 
nents have their origin in the interstellar medium surround¬ 
ing the nebula and therefore, we can assume that these fore¬ 
ground fluxes are similar on the whole nebula NGG 6888. 
By comparing the total fluxes from spectrum A with these 
values it is found that the foreground contribution is only 
around 1 per cent for all the identified lines, except for the 
faint [O IIlJAA 4959,5007 lines where this flux contribution 
reaches 5 per cent. We can therefore assume that the inte¬ 
grated fluxes of the spectrum at this position represent well 
the nebular behaviour. 

From the data in Table |2] we determined the electronic 
temperature Te and density rie from the available diagnos¬ 
tic ratios. In the case of positions B and C only the intense 
blueshifted components were used. In all the cases, only Te 
derived from the [N II] line flux ratio could be determined 
because the [O III] A4363 line (crutial for determining Te [O 
III]) was not detected in the spectrum. The rie value comes 
from the usual [S II] diagnostic ratio because the lines com- 
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Table 4. Physical parameters and ionic and total abundances 


Parameter 

Position A 

Position Bf, 

Position Cfe 

EV92 Position 

Position B2 (FM12) 

M-D14 

Te[0 III](K) 




9500: 

(6855+639|^ 

9050+810 

Te[N II](K) 

78001^“ 

11700+3000 

9100+1150 

8400+100 

8537+527 


UefS II] (cm“3) 

240tlll 

1801 ? 75 ° 

91 c;+580 
^-*-0-130 

371 

108+26 

480: 

UefO II] {cm-3) 

-*-^^-110 



301 


310+30 

12+log 0+/H+ 

o -I q + 0.37 



7.54+0.34 

7.96+0.14 

7.45+0.14 

12+log 0++/H+ 

7 7Q+0.28 
' • '^-0.16 

7.80t°-j6 

7 qq+0.26 

7.98+0.25 

8.57+0.19 

8.02+0.11 

12+log N+/H+ 

o Q -| + 0.21 

7 -| q + 0.53 
' •-*^^- 0.21 

7 7 Q + 0.20 

'•'^-0.14 

7.83+0.30 

7.63+0.08 

7.54+0.08 

12+log S+/H+ 

6 34+'^-21 

3 Qf;+0.51 

^•'^'^-0.19 

5 66+0-21 

^•^^-0.13 

5.71+0.23 

5.80+0.08 

5.43+0.09 

12+log S++/H+ 


6 45+0-50 
0.40_q 5Q 

6 75+0-35 
'^-0.40 




12+log Ar++/H+ 

'^•^^- 0.12 

O.UU_Q 20 

'^•-*-^-0.13 



6.36+0.10 

12+log He+/H+ 

11.20+0.18 

11.29+0.14 

11.21+0.36 

11.24+0.06 

11.20+0.06 

11.21+0.03 

12+log O/H 

8 - 29 ^° u 9 



8.11+0.28 

8.66+0.16 

8.20+0.09 

12+log N/H 




8.40+0.35 

8.34+0.23 

8.54+0.20 

12+log S/H 

7 -| c + 0.20 
' --*-^- 0.12 





6.77+0.20 

12+log Ar/H 

5 47 + 0.20 
'^•^'- 0.12 





6.41+0.11 

log N/O 

0 19+0-^^ 
^•-^^-0.32 



0.29+0.30 

-0.33+0.28 

0.22+0.22 


Temperature derived from 


Perez-Montero & Contini 


(2009 I formula. 


Table 5. Atomic data 



Transition 

Collisional 

Ion 

probabilities 

strengths 


N+ 

0 + 

0 ++ 

s+ 


Galavis, Mendoza &: Zeipper 

119971 

Wiese, h’lihr Deters 

(pm 

Zeippen (|l982 ' 


Wiese, i-'uhr ^ Deters 

(1996 


Wiese, i^'uhr ^ Deters 

(1996 


ytorev ^ Zeippen (|2U0U ' 


Tayal I 2011 1 


Pradhan et al. 

(2006 


'i'aval 




Aggai 

wal Keenan 

1999 


S++ 

Podobedova, Kelleher Wise 

;ioo9 


Galavis, Mendoza ^ Zeippen 

1995 


Ar++ 

Kautman bugar 

(p^H6| 


Galavis, Mendoza Zeippen 

199S 

1 


Tayal fc Zatsarinny i 2010 1 


ing from this ion are present in all the spectra. The tempera¬ 
ture and density as well as the ionic abundances, derived by 
using these physical conditions and the line intensities from 
Table 2, are shown in Table All these quantities as well as 
the errors were determined with the software PyNeb 1.0.Ib9 


(Luridiana, Morisset & Shaw 20151 using the atomic data 
shown in Table The total abundances were only derived 
for position A where [O II] AA3726,3729 were detected. The 


ionization correction factors (IGF) by Kingsburgh & Barlow 
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( 1994|) were used with this purpose. The results are listed in 
Table 0 

In addition to our results, in Table 0 we are including 
the results obtained by EV92 for their blueshifted compo¬ 
nent, and the data by FM12 corresponding to the inner zone 
observed by them, i.e., that labeled in their paper as B2. 
Both zones (EV92, and B2 EM12) are clearly not associ¬ 
ated with the rim region of NGC 6888 as it is shown in our 
Fig. 0 where their positions are marked. Also the results by 
M-D14, for their position in the rim, have been listed. 

From these data we can make a direct comparison of 
the chemistry in different zones of NGC 6888, as obtained 
by different authors. 

In Table 0 it is found that our values of Te [N II] elec¬ 
tronic temperature and density for positions A, B, and C 
agree with the values derived by EV92, FM12 and M-D14. 
According to our observations, the plasma in position B(, has 
the highest electronic temperature (although with a large 
uncertainty); similar values were reported by FM12 for the 
gas located farther than the edge of NGC 6888 and although 
their regions and our B position are not spatially close an 
explanation for the high temperature in position B could be 
that the main nebular emission comes from a low-density 
component of NGC 6888 if we assume this ring nebula is 
made of two gaseous components (see §4.2.2). 

The abundances we obtain for position A in the rim 
are in good agreement with the previously published abun¬ 
dances (except for position B2 of FM12, see below). The 
fact that the values of O/H and N/H are compatible within 
the error bars for spatially different regions of the nebula 
indicates that the hypothesis of a chemically homogeneous 
nebulae need to be explored. 

The abundances found for regions slightly farther than 
the edge of NGC 6888 (positions Bl, El, E2 and MBS of 
FM12) are not taken into account because, according to 
FM12, there is no nitrogen enrichment there as this com¬ 
ponent was probably ejected during an early evolutionary 
stage of the star WR 136 as proposed by |Garcfa-Segura,| 


Langer & Mac Low (19961. 


The region B2 of FM12 presents an O/H abundance 
much larger than the values of the other positions. Note 
that FM12 did not obtain directly the electron temperature 
for the [O III] zone in their B2 region, but rather defined it 
using the relation (2) from Perez-Montero & Contini (20091 
which was derived for other kind of objects as galaxies and 
H II regions, and given that NGC 6888 is an extented neb¬ 
ula which has a complex structure, probably made of two 
gaseous components, the changes in temperature we detect 
across the nebula cannot be assumed to follow the same be¬ 
havior as found for global observations, since in addition we 
must deal with the ionization structure. 

Therefore, their very low [O III] temperature value (not 
found by EV92 who directly determined it in a region close 
to B2) explains the high value they consequently obtain for 
ionic abundance and their high total O/H. If we 
use the Te[0 III] of EV92 the abundance is only 8.01 
and the total abundance of oxygen is 8.29 which is in good 
agreement with the rest of the values shown in Table 0 

In addition FM12 obtain a low log N/O ratio of —0.33 
for this zone, which is also questionable, as this is contrary 
to the chemical enrichment predicted (log N/O ~0.20) by 
the stellar evolution models (discussed below), where N- 



Figure 5. Scheme of the geometry assumed to build the pho- 
toinization model of NGC 6888. We can see the high-density 
and clumpy component (Fmodg) as an ellipsoide covered by 
the spherical low-density component [Fmod]^). In the model, the 
edges of the ellipsoide match the size of the observed semi-axis of 
NGC 6888. 


enrichment should be detected also near the star. This is 
actually the case for the also inner position observed by 
EV92 who obtained log N/O = 0.29. Note that the N/O 
abundance ratio is defined as equal to the N"'"/0^ ratio 
(ICF=1), while we will see in §5.3 and Fig. 0 that this is 
not the case for observations obtained through small slits in 
the inner zone of this extended nebula. 

The values by FM12 could be indicating a chemically 
inhomogeneous nebula, as these authors suggested. However 
given the above discussion, it is possible instead to assume 
that the main structure of NGC 6888 has a homogeneous 
chemical abundance. The existence of such homogeneous 
nebulae are supported by the stellar evolution models for 
WR stars. For example, from Figure 16 of |Toala &: Arthur| 
( |2011[ ) we can see that along the evolution of a WR star, the 
atmospheric log N/O ratio reaches a high value of 0.20 in 10 
000 years and remains almost constant for around 300 000 
years, in the case of a non-rotational star with initial mass 
of 40 Mq. 

In a similar way M-D14 followed the time evolution of 
the N/O ratio for NGC 6888, concluding that the material 
of the rim was ejected during the Red Super Giant phase 
highly nitrogen-enriched and the subsequent WR phase did 
not provide additional enrichment (see their Figure 5). 

Therefore, given that the dynamical age of NGC 6888 


has been estimated to be of about 20,000 - 40,000 yr (van 


Marie, Langer & Garcfa-Segura 20051, the idea of a chem¬ 


ically homogenous nebula for NGC 6888 is consistent with 
the evolutionary models of a massive star becoming a WR 
star. 

In the following we explore this possibility, by means of 
a photoionization model. 
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Table 6. Abundance from the nebular and stellar models in 
12+log(X/H) 


Element 

Nebular 

starS 

He 

11.21 

12.25 

C 

8.86 

7.84 

N 

8.40 

9.95 

O 

8.20 

8.76 

Ne 

7.51 


Si 


8.63 

P 


8.01 

S 

7.10 

6.34 

Ar 

6.41 


Fe 

6.33 

8.35 


4-2.2 Photoionization model of NGC 6888 


To reproduce the optical spectra of NGC 6888 we computed 
a photoionization model using the pyCloudy package (Moris- 
set 2013) based on the ID photoionization code Cloudy (Fer- 
land 2008). It allows us to generate 2D images for each emis¬ 
sion line from a pseudo-3D model and to apply a mask on 
them, reproducing the effect of the Hnite aperture through 
which the observations were performed. 

For our nebular model we used the version CIO.O of 


Cloudy, last described by Ferland et al. | 2013[). We assumed, 
as it is apparent from the optical images | |Gruendl et al. 
2000 Moore, Hester fc Scowen||2000 1, that the NGC 6888 
projected morphology is the result of two gaseous compo¬ 
nents, as shown in Figure one bright, granular region 
emitting mainly low ionization lines (e.g., [N II]), and the 
other one, more diffuse, mainly emitting the higher ionized 
emission lines (e.g., [O III]). 

The denser (brighter) component has an ellipsoidal 
shape with an inner minor semi-axis of 80 arcmin, a mi¬ 
nor to major semi-axis ratio of 0.33 with a H density of 400 
cm“®; the other component has a lower density of 1 cm“^, 
and a spherical geometry with a central cavity with a radius 
of 20 arcmin. In addition we set a constant density prohle 
for each component, with a filling factor /mod^, = 1.0 for the 
low-density component and fmod„ = 0.1 for the high-density 
component. 

As we aim to build a consistent model of the system 
we used as ionizing sources the stellar atmosphere models 
of the star WR 136, as calculated with the CMFGEN code 
and described in the previous section. We therefore ensure 
to have well constrained stellar parameters, the ionization 
parameter and the ionization degree of the nebula through 
the stellar luminosity and temperature. 

In a simple way, we calculate a chemically homogeneous 
model for NGC 6888 and we verify if observations can be 
reproduced. 

The chemical abundances were hrst set to the values de¬ 
rived for position A in §4.2.1 (Table|^ which were completed 
with the abundances of C, Ne and Fe published by M-D14. 
Then the abundances were slightly modified until the pho- 
toionized model reproduces the observed line intensities.The 
hnal values which describe better the chemical composition 
in NGC 6888 are given in Table 

The coordinates for the slit positions, from where the 
model integrated fluxes were obtained, were set in angular 


units with respect to the central star adopting a distance of 
1.45 kpc for NGC 6888 ( van Leeuwen|[2007 1. In the optical 
images the star does not seem to be in the geometrical center 
of the nebula, but as we are only modelling the northwest 
part of it, the geometry of the model is an ellipsoid with the 
geometrical center on the central star position and whose 
edge matches the observed region. 

The final intensities of the emission lines are obtained by 
summing up the contributions of each component, assuming 
that the high-density region weight on the total intensity is 
w. Thus: 


Pmod — fmodj^Pmodj^ll^ rc) “t“ fmodFmod^'O) 


(3) 


where F-modL Fmodn the line intensity values from 
the low-density and high-density components, respectively. 
Once the abundances were fixed, only this weight value was 
changed, from one position to another, to fit the observed 
fluxes. 

The comparison between the model predictions and 
the observed dereddened line fluxes is presented in Table 
and Table corresponding to the photionization model 
results by using star! and star3 SED’s as ionizing source, 
respectively. The q parameter, in these tables, shows the 
goodness of our model and is defined following [Morisset fc| 
Georgiev (20091 as: 


loglF-fjioR log^Foifs') 

t 


(4) 


where F^od is the value returned by the model. Fobs is the 
observed flux and t is the accepted tolerance in dex of this 
observable. We define the tolerance to be 


' =(‘ + xi?) 

where SFobs is the measured error. When a line flux is well 
reproduced the q value is between —1.0 and 1.0. 

The results in Tablej^show that, under the assumptions 
made for this model and using star! as ionizing source, the 
observed line fluxes are not well reproduced. In particular 
the lines from highly ionized elements (O^"*", Ar’’'^) are ex¬ 
tremely overpredicted. The best case is for the values in 
position A where only half of the whole line fluxes (the low 
excitation ones) detected in this region are well reproduced. 
The model results are worse for the other positions. 

The fact that the ionization of the nebula is strongly 
overpredicted could be due to: 1) a too high ionization pa¬ 
rameter 2) a too high stellar effective temperature. No¬ 
tice that this model was obtained using only the high den¬ 
sity component (w = 1). Adding any contribution of the low 
density component would made the situation even worse. 

As R, Qs and uh are well constrained through the HIP- 
PARCOS distance to the object, its observed angular size, 
the observed stellar luminosity and the nebular density com¬ 
ing from the plasma diagnostic, then it follows that U is also 
well constrained. The only remaining way to reproduce the 
observed low value of the [O IIl]/]0 II] line ratio is to con¬ 
sider a star with a temperature lower than the value derived 


1/ = 4 ^^* where Q* is the ionizing photon rate, ng the hy¬ 
drogen number density and R the distance to the ionizing source. 
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Table 7. Comparison between the deredden observational and model fluxes (normalized to = 100) using starl as ionizing source 

and w = 1 


A (A) 

Ion 

Fobsyl 

FmodA 

OA 

Fobss 

Fmods 

qt 

Fobs(7 

Fmodc 

qc 

3726.03 

[O II] 

52.6±14.3 

52.5 

0.00 







3728.81 

[O II] 

62.6±9.5 

61.7 

-0.13 







3970.07 


16.8±6.1 

15.1 

-0.35 







4340.46 


47.1±5.7 

46.3 

-0.13 

48.8±13.4 

46.8 

-0.17 

47.7±3.5 

47.0 

-0.17 

4958.91 

[O III] 

23.8±2.4 

80.6 

14.9 

97.0±6.2 

121.5 

2.91 

52.1±4.5 

136.0 

12.51 

5006.84 

[O III] 

67.9±5.6 

243.7 

18.2 

289.2±16.3 

365.6 

3.51 

154.2±9.2 

409.2 

18.81 

5015.67 

He I 

4.7±1.8 

2.8 

-1.74 

6.0±1.8 

4.2 

-1.62 

4.7±3.1 

4.2 

-0.21 

5754.64 

[N II] 

3.8±1.2 

4.5 

0.74 

2.3±1.3 

1.3 

-1.33 

2.3±0.9 

1.1 

-1.22 

5875.61 

He I 

19.9±2.3 

22.4 

0.96 

21.3±2.5 

23.22 

0.68 

23.5±2.1 

23.2 

-0.14 

6312.10 

[S III] 




2.2±1.3 

4.4 

1.67 




6548.04 

[N H] 

168.4±11.5 

125.8 

-4.25 

37.1±2.8 

32.1 

-2.83 

70.9±3.5 

24.3 

-15.70 

6562.80 

H<, 

285.9±18.3 

291.4 

0.30 

286±14.1 

290.1 

0.16 

286±10.8 

289.2 

2.00 

6583.46 

[N H] 

502.1±31.9 

371.1 

-5.10 

106.8±6.1 

94.7 

-2.59 

227.0±9.7 

71.8 

-19.13 

6678.15 

He I 

6.4±1.2 

6.3 

-0.13 

6.4±2.8 

6.6 

0.10 

5.2±0.7 

6.5 

2.79 

6716.44 

[S II] 

22.9±2.5 

31.3 

3.92 

7.0±2.1 

10.1 

1.66 

7.4±0.8 

8.0 

0.90 

6730.81 

[S II] 

19.9±1.9 

28.4 

4.73 

5.7±2.4 

9.8 

1.73 

6.3±1.2 

7.8 

2.36 

7065.17 

He I 

3.2±1.0 

3.9 

0.84 

6.0±2.6 

4.4 

-0.88 

6.4±0.8 

7.83 

0.66 

7135.80 

[Ar HI] 

9.2±1.1 

25.4 

5.56 

15.6±2.7 

25.5 

0.52 

lO.Oil.l 

26.3 

5.64 


Table 8. Same as Tab.[^ but using star3 as ionizing source 


A (A) 

Ion 

Fobs^ 

Fmod^ 

qA 

Fobsg 

Fmodg 

1b 

FobS(;7 

Fmod(;7 

qc 

3726.03 

[O H] 

52.6±14.3 

50.6 

-0.13 







3728.81 

[O H] 

62.6±9.5 

60.0 

-0.37 







3970.07 

H, 

16.8±6.1 

15.1 

-0.36 







4340.46 

H.^ 

47.1±5.7 

45.7 

-0.27 

48.8±13.4 

43.8 

-0.45 

47.7±3.5 

45.1 

-0.82 

4958.91 

[O HI] 

23.8±2.4 

22.04 

-0.81 

97.0±6.2 

68.0 

-3.03 

52.1±4.5 

51.5 

-0.13 

5006.84 

[O HI] 

67.9±5.6 

66.33 

-0.30 

289.2±16.3 

204.4 

-3.45 

154.2±9.2 

155.0 

0.12 

5015.67 

He I 

4.7±1.8 

2.74 

-1.70 

6.0±1.8 

1.8 

-4.75 

4.7±3.1 

2.9 

-0.93 

5754.64 

[N H] 

3.8±1.2 

3.63 

-0.12 

2.30±1.3 

1.0 

-1.84 

2.3±0.9 

2.0 

-0.19 

5875.61 

He I 

19.9±2.3 

22.7 

1.21 

21.3±2.5 

19.8 

-0.62 

23.5±2.1 

21.8 

-0.97 

6312.10 

[S HI] 




2.2±1.3 

1.4 

-0.81 




6548.04 

[N H] 

168.4±11.5 

163.2 

-0.43 

37.1±2.8 

37.1 

0.00 

70.9±3.5 

71.7 

0.14 

6562.80 

He 

286±18.3 

291.7 

0.29 

286±14.1 

296.2 

0.37 

286±10.8 

293 

0.35 

6583.46 

[N H] 

502.1±31.9 

481.7 

-0.71 

106.7±6.1 

109.6 

0.70 

227.0±9.7 

211.48 

-1.10 

6678.15 

He I 

6.4±1.2 

6.44 

0.04 

6.4±2.8 

5.6 

-0.33 

5.2±0.7 

6.2 

2.13 

6716.44 

[S H] 

22.9±2.5 

15.87 

-3.51 

7.0±2.1 

2.1 

-4.58 

7.4±0.8 

3.9 

-6.75 

6730.81 

[S H] 

19.9±1.9 

14.6 

-3.28 

5.7±2.4 

2.0 

-2.85 

6.3±1.2 

3.8 

-4.86 

7065.17 

He I 

3.2±1.0 

3.6 

0.41 

6.0±2.6 

3.3 

-1.78 

6.4±0.8 

3.7 

-1.64 

7135.80 

[Ar HI] 

9.2±1.1 

10.5 

1.14 

15.6±2.7 

8.79 

-0.15 

lO.Oil.l 

10.2 

0.24 

w (see text) 

7.6 X 10-2 



1.7 X 10“"^ 



3.8 X 10-"* 




for model starl, i.e., we must not deal with the number of 
ionization photons but rather with their energy. 


The fit improves very much when we use the model 
star3 which has T, = 70 000 K. The ionization of the dense 
gas is even lower than observed, leading to the full use of our 
hypothesis of the 2-components model. The predicted emis¬ 
sion lines from the corresponding photoionization model are 
compared with observations in Table It can be seen that 
almost all the fluxes are well reproduced, indicating that the 
star WR 136 emits the UV ionizing photons corresponding 
to a star with this low effective temperature. 


5 DISCUSSION 

5.1 Effective temperature and mass-loss rate of 
WR 136 

As it can be seen in Figure]^ and it was discussed in §4.1, 
at least three stellar models, {starl, star2 and star3) with 
different T, and therefore a different Rt fit similarly well the 
stellar observations, in particular the stellar continuum and 
several emission lines. The differences among these models 
appear mainly in the He lines, where the He Il/He I line 
ratio is well reproduced when considering an effective tem¬ 
perature of 110 000 K but it is understimated when the 
effective temperature is reduced, i.e., at a lower stellar tem¬ 
perature the stellar He I lines are stronger and the stellar 
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He II lines become weaker than observed. As we discussed 
before, it seems that there is not a unique solution for the 
stellar model. 

In Figure we show a comparison between the UV re¬ 
gion of the stellar models from Table We also overplot the 
UV behaviour of two blackbodies with temperatures of 70 
000 and 46 000 K which correspond to the temperatures T* 
and r 2/3 of model starS. 

We can see that the key difference between the stellar 
models comes from the far UV region of the spectra, whose 
UV photons are absorbed by the nebula (between 1.0 and 
~ 7.0 Ryd). This difference is well appreciated when we 
compare the models with the blackbodies. From this figure 
it is clear that a blackbody with the same temperature than 
the modelled spectra, 70 000 or even 46 000 K, shows an 
excess of UV photons which would lead to a much more 
ionized nebula. It is also evident the high absorption of the 
UV photons emitted by the star, due to the presence of the 
stellar wind. 

Owing to there are not stellar observations in the far 
UV to discriminate which stellar model is reproducing the 
real spectral energy distribution of the star WR 136, it is 
necessary to add an external restriction. In this case, such 
additional restriction is the associated nebula NGC 6888 
that is highly sensitive to this UV ionizing radiation. 

In addition, the combination of the temperature and 
mass-loss effects in the stellar wind, considered for this 
model, can reproduce the line profiles coming from recom¬ 
bination process. Therefore, if we are close to a good value 
for T*, it is the same for the mass-loss rate, M. This is 
important because M represents a good constrain for the 
hydrodynamical simulations. 

As in this study we have found a well constrained value 
for M (log M = —4.95), which is similar to the values found 
by other authors (Table [3|, we c an rule out the assumption 
made recently by |Zhekov fc Parl^ ( |2011[ ) where they conclude 
that it is necessary a M 3 - 4 times lower than observed for 
this type of stars to reproduce the X-ray flux of the emitting 
hot bubble. 


5.2 Is NGC 6888 mainly photoionized? 


The presence of a very fast stellar wind coming from the star 
WR 136 could imply that part of the ionization of the gas 
is due to shocks. Diagrams made by|Veilleux & Osterbrock| 
|T987| ) and |Baldwin, Phillips fc Terlevich| ( |1981[ ) are usually 
employed to determine if a line emitting region is mainly 
ionized by OB stars (like H II regions) or by other source 
(AGN, LINERS, shocks). See for example |Allen et al.| ( |200^ 
for the case of shocks. 

In Figure we show an [OIII] 5007/H/3 vs. [Nil] 
6584/Ha diagram. The solid line shows the separation be¬ 
tween global photoionization models for Star Forming re¬ 


gions (NSF galaxies) and AGNs as defined by Kauffmann et 
[aL] ( (20031 ). The dashed line represents the same separation, 
but shifted on the log (]N II] 6584/Ha) axis by 0.96 dex. 


which is the difference between the N/O adopted by Kauff- 
(mann et al.| ( |2003| ) and the N/O of NGC 6888 deduced in this 
work. The three black diamonds correspond to observations 
for positions A, B and C and are obviously located to the 
right side of the solid line, but well inside the photoionized 
zone with N/O ratio similar to NGC 6888 value. 



Figure 6. Comparison of the ionizing region of the spectra cor¬ 
responding to stellar models with diferent effective temperature. 
The flux is not dereddened and is shown as emitted at a distance 
of a 1 kpc. Colors are as in Figure|^ The black solid line and the 
dashed line correspond to blackbodies of 70 000 and 46 000 K, 
respectively, according to the temperatures T, and T 2/3 of model 
starS. 


The effect of a limited aperture size, which is not in¬ 
cluding the whole nebula is illustrated by the position of the 
color points in our Fig.[^ They show the emission line ratios 
of every spaxel of our 2D images obtained from our models. 
The colors correspond to different values of w (the weight 
of the dense component), the blue symbols are obtained for 
the case of considering only the low-density component, and 
the red ones for the case of considering only the high-density 
component, the green symbols being an intermediate situa¬ 
tion. As we can see, the points spread in an ample region of 
the diagram. The black square symbol is the position in the 
diagram of the global model, as if it were fully included in 
the slit. 

Therefore if just a small portion of a whole nebula is 
observed, one must be careful in using diagnostic diagrams 
to decide if the emission of a nebula is caused by the pres¬ 
ence of shocks, because as we demonstrated in Fig. 0 only 
photoionization and an enhanced N abundance are sufficient 
to obtain points located out of the regions represented by 
global photoionized models (solid line). 


5.3 Oxygen and nitrogen abundances in NGC 
6888 

To deeper investigate our assumption of chemical homogene¬ 
ity of NGC 6888, in Fig. we added the observations by 
EV92, M-D14. It is found that they lie inside or close to 
our model points. Also we added the observations by FM12 
(marked as black dots and stars) resulting from their ID 
analysis of different regions of NGC 6888. We can see that 
their points show a large spread in this diagram. To explain 
the behaviour of their observations and to derive values of 
the N and O abundances, FM12 performed a grid of pho¬ 
toionization models using the code Cloudy varying the ion¬ 
ization parameter and the N/O ratio. As they used global 
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Figure 7. Plasma diagnostic diagram [O IIIJ/H/3 vs. [N II]/Ha. The solid curve corresponds to the diagram in |KaufFmann et al.| l [200^ 
and the dashed curve is the same one but shifted by 0.96 dex given the difference between the N/O values in that work and in ours. For 
three different w values, all the points that make the photoionization model of NGC 6888, using star3 as ionization source, are shown. 
The observed ratios are marked as black diamonds and the black square represents the values of integrated fluxes coming from the whole 
nebula. The observations by FM12 are also shown: stars represent the N-poor regions in the immediate outer zone of NGC 6888 while 
the black points are observations of the main structure of the nebula with N-rich abundance. Observations by EV92 and M-D14 are also 
show in the diagram as black triangles. 


models for each observation without taking into account that 
they were observing only a part of the nebula they concluded 
that there are important differences in the N/O ratio across 
NGC 6888. 

However, our figure shows that this is not necessarily 
the case. In particular, let us focus on the data of FM12 
marked as black dots, which correspond to their regions in¬ 
side NGC 6888 with large nitrogen enrichment. We can see 
that all these points, but the point labeled X2 in their pa¬ 
per, are close to our model points and are following their 
trend. The rest of the points (stars) show a large dispersion 
compared to our model points and it is not expected for the 
model to reproduce their positions because they belong to 
regions located out from the edge of NGC 6888 with low 
(solar) nitrogen abundance. 

Therefore, Figure illustrates the fact that all the ob¬ 
served points within the nebula can be reproduced by a 
rather simple chemically homogeneous model. 

To continue the discussion about the chemical homo¬ 
geneity in NGC 6888, we analyse the behaviour of the IGF 
used to determine the N/O ratio. Let us remember that it 
is generally assumed that N"''/0^ is equal to N/O. 

From our photoionization model, we are able to deter¬ 
mine in each region of the nebula the IGF needed for the 
N"''/0’*’ ratio to obtain the N/O value. We show in Figure!^ 
a map of this IGF for three different values of w. 10“®, 10“ 
and 10“^. It is evident that for a given w, this IGF varies 


over the nebula showing a difference as large as a factor 
of two. In adittion, we can see that the value of this IGF 
changes as the contribution of the high-density component 
is modified. 

For example, in the case when w = 10“^, the IGF for 
position A is around 1.0. In general, we can see in Figure 
that for positions closer to the star the IGF is around 1.5 
or larger. Therefore, in the zones near to the star an IGF 
greater than one is needed to estimate the N/O ratio. Hence, 
this indicates that N/O could be understimated when it is 
simply assumed that N"'"/0^ is equal to N/O. In this way, 
to calculate the nitrogen abundance for position B2 of FM12 
and for EV92 a better value for this IGF is 1.5. Once this 
IGF is used for these positions the N/O ratio increases rel¬ 
ative to the values in Table to 0.47 for EV92 and —0.15 
for B2 FM12. Therefore these values are still discrepant. In 
particular, for the position B2 the N/O ratio reaches a value 
of —0.15 lower than the others for the rim, which is contrary 
to the predictions of stellar evolution models. 


The differences seen of the IGE in different zones of the 
map are due to the combination of two effects: the ionization 
structure of the nebula and the relative contribution to the 
fluxes of the high density clumpy region. 


A similar problem with the IGFs is affecting the results 
of Stock & Barlow ( 2014[ ) who used Herschel observations 
of the [O III] 88 and [N III] 57 fj,m lines, emitted in two 
positions of NGG 6888, to determine the N/O abundance 
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Figure 8. Map of the ICF for the N+ZO"*" ratio from a quadrant 
of our best photoionization model which is chemically homoge¬ 
neous. The star is assumed to be at the origin of the map. Three 
different values of w were used to show how the values of the ICF 
changes over the whole structure of nebula. Our positions and the 
position of FM12 are shown. 


ratio under the assumption that is equal to N/O. 

For NGC 6888 this is a very rough approximation because 
only a small fraction of the nebula was included in their 
slits and Stock & Barlow (20141 did not take into account 
the nebular ionization structure, which indicates that the 
ion N^"*" is only dominant over in the zones near to the 
star. 


The previous discussion illustrates very well that de¬ 
tailed modelling (even the simplest as we do here: 2 almost 
roundish homogeneous components) is a very powerful way 
to identify effects that, if not considered, can easily drive 
the observer to wrong conclusions. 


6 CONCLUSIONS 

We computed several atmosphere models using the code CM- 
FGEN to reproduce the stellar observations of WR 136. The 


characteristics of the star such as L, M and 14o agree with 
the values found by other authors, i.e., HGL06 and|Crowther| 


& Smith (19961. We also derived the stellar abundances of 
He, G , N, O, Si, P, S and Fe (listed in Table|^ which agree 
with HGL06’s model. 


We find that models with very different effective tem¬ 
perature, Tt from 70 000 K to 110 000 K, can reproduce 
fairly well the stellar continuum and several emission lines, 
although the He lines, indicative of temperature, are worse 
reproduced by the cooler models. Therefore there is a de¬ 
generacy in atmosphere models, than apparently, as pointed 
out by |Hillier| ( |1991| , is due to our misunderstanding of the 
structure of the wind velocity field. 

To decide which stellar model produces the correct UV 
amount of ionizing photons coming from WR 136, we com¬ 
puted a photoionization model for the associated ring nebula 
NGC 6888, using the SED produced by the stellar models. 
This led us to conclude that WR 136 follows better the spec¬ 
tral behaviour of a star with a T* of 70 kK. 


The photoionization model was built assuming chemical 
homogeneity and including two components: a low-density 
region enclosing a denser region. 


When these two gaseous components are considered and 
assuming only photoionization, we were capable to repro¬ 
duce the observed ionization structure in both shells, con¬ 
trary to the presence of shocks as the probable mechanism 
producing the observed ionization structure of NGC 6888. 
We cannot immediatly rule out the importance of shock ef¬ 
fects, but given that some diagnostic diagram such as the 
[O III] 5007/H/3 vs. [N II] 6584/Ha indicate emission of H II 
type, and that we can reproduce several emission features of 
the nebula by only assuming photoionization, we established 
that the photoionization heating is more important than the 
presence of shocks. This latter idea was already mentionned 
by EV92 combining the diagnostic diagrams of |Sabbadin,| 
[Minello fc Bianchini| ( |1977| ) with their observations in a zone 
close to the central star and shown in Figure]^ 

The nebular abundances derived from our photoioniza¬ 
tion model are similar to the values derived from our position 
A and the values reported by EV92 and M-D14. 

It is worth to notice that both the stellar atmosphere 
of WR 136 and the nebular plasma present a large N- and 
C-enhancement and an O depletion, as it can be seen in 
Table This contributes to the idea that the observed gas 
was ejected by the central star through the mass-loss driven 
by the stellar wind. 


Considering only photoionization and a chemically ho¬ 
mogeneous nebula, we were able to explain the pretended 
differences found in the N/O ratio (derived from N"'"/0^ or 
N++/0++ in different zones of the nebula) by FM12 and 
[Stock fc Barlow (20141, who proposed a chemically inho¬ 
mogeneous nebula. These pretended differences are due to 
effects of the ionization structure within the nebula. With 
our simple model we show that the ICFs for global nebulae 
should not be use for observations of small zones observed 
in extended nebulae. 


From our results we have learned that there are sev¬ 
eral stellar models with different set of values that could 
reproduce the observables, like the models starl and star3 
presented in this study. Therefore external constraints are 
crucial to determine which atmosphere models represent the 
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real stellar properties. The nebulae around stars can help to 
deal with the degeneracy mentioned above. 

Finally, in this work we show the advantages of mod¬ 
elling simultaneously the stellar and nebular spectra of a 
Wolf-Rayet star plus nebula system. We verify whether a 
stellar atmosphere model of the star WR 136 provides a 
good description of the stellar observations and is able to 
satisfactorily reproduce the nebular characteristics such as 
the ionization degree. This allow us to better constraint all 
the physical parameters of both the star and the nebula, 
including their physical conditions and their chemical abun¬ 
dances. 
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